numerical integration of x^3/(e^x-1) from 0 to pi